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Abstract 

The novel Parton-Hadron-String Dynamics (PHSD) transport approach is applied 
to nucleus-nucleus collisions at RHIC energies with respect to differential hadronic 
spectra in comparison to available data. The PHSD approach is based on a dynam- 
ical quasiparticle model for partons (DQPM) matched to reproduce recent lattice- 
QCD results from the Wuppertal-Budapest group in thermodynamic equilibrium. 
The transition from partonic to hadronic degrees of freedom is described by covari- 
ant transition rates for the fusion of quark-antiquark pairs or three quarks (anti- 
quarks), respectively, obeying flavor current-conservation, color neutrality as well as 
energy-momentum conservation. Our dynamical studies for heavy-ion collisions at 
relativistic collider energies are compared to earlier results from the Hadron-String 
Dynamics (HSD) approach - incorporating no explicit dynamical partonic phase - as 
well as to experimental data from the STAR, PHENIX, BRAHMS and PHOBOS 
collaborations for Au+Au collisions at the top RHIC energy of \fs = 200 GeV. 
We find a reasonable reproduction of hadron rapidity distributions and transverse 
mass spectra and also a fair description of the elliptic flow of charged hadrons as 
a function of the centrality of the reaction and the transverse momentum px- Fur- 
thermore, an approximate quark-number scaling of the elliptic flow V2 of hadrons is 
observed in the PHSD results, too. 
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1 Introduction 



Present experiments at the Relativistic Heavy-Ion Collider (RHIC) have reached for 
short time scales the conditions met in the first micro-seconds in the evolution of the 
universe after the 'Big Bang'. The 'Big Bang' scenario implies that on these time 
scales the entire state has emerged from a partonic system of quarks, antiquarks and 
gluons - a quark-gluon plasma (QGP) - to color neutral hadronic matter consisting of 
interacting hadronic states (and resonances) in which the partonic degrees of freedom 
are confined. The nature of confinement and the dynamics of this phase transition 
is still an outstanding question of todays physics. Early concepts of the QGP were 
guided by the idea of a weakly interacting system of massless partons which might 
be described by perturbative QCD (pQCD). However, experimental observations at 
RHIC indicated that the new medium created in ultrarelativistic Au+Au collisions 
is interacting more strongly than hadronic matter (cf. [1] and Refs. therein). It is 
presently widely accepted that this medium is an almost perfect liquid of partons 
[2-7] as extracted experimentally from the strong radial expansion and the scaling of 
the elliptic flow v 2 (pt) of mesons and baryons with the number of constituent quarks 
and antiquarks [6,8,9]. In order to explore the transport properties of this partonic 
medium microscopic studies based on non-equilibrium dynamics are mandatory. 

A consistent dynamical approach for the description of strongly interacting systems 
- also out-of equilibrium - can be formulated on the basis of Kadanoff-Baym (KB) 
equations [10,11] or off-shell transport equations in phase-space representation, re- 
spectively [11-14]. In the KB theory the field quanta are described in terms of dressed 
propagators with complex selfenergies. Whereas the real part of the selfenergies can 
be related to mean- field potentials (of Lorentz scalar, vector or tensor type), the 
imaginary parts provide information about the lifetime and/or reaction rates of time- 
like 'particles' [4] . Once the proper (complex) selfenergies of the degrees of freedom 
are known, the time evolution of the system is fully governed by off-shell transport 
equations (as described in Refs. [11-14]). The determination/extraction of complex 
selfenergies for the partonic degrees of freedom has been performed before in Refs. 
[4,15,16] by fitting lattice QCD (1QCD) 'data' within the Dynamical QuasiParticle 
Model (DQPM). In fact, the DQPM allows for a simple and transparent interpre- 
tation of lattice QCD results for thermodynamic quantities as well as correlators 
and leads to effective strongly interacting partonic quasiparticles with broad spec- 
tral functions. We stress that mean-field potentials for the 'quarks' and 'gluons' as 
well as effective interactions can be extracted from 1QCD within the DQPM (cf. 
Ref. [16]). For a review on off-shell transport theory and results from the DQPM in 
comparison to 1QCD we refer the reader to Ref. [14]. 

In preceding works [17,18] two of the authors have presented first Parton-Hadron- 
String-Dynamics (PHSD) transport calculations for expanding partonic fireballs as 
well as nucleus-nucleus collisions at Super-Proton-Synchrotron (SPS) energies of 40 
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to 160 A GeV. The studies in Ref. [17] have addressed expanding partonic fireballs 
of ellipsoidal shape in coordinate space that hadronize according to local covariant 
transition rates. It was found that the resulting hadronic particle ratios turn out 
to be in line with those from a grandcanonical partition function at temperature 
T PS 170 MeV rather independent from the initial temperature of the partonic system. 
Furthermore, the scaling of elliptic flow with initial spatial eccentricity indicated a 
dynamical evolution of the system close to ideal hydrodynamics, which so far has 
been successfully employed for the description of experimental data at RHIC [19- 
22] 1 . Additionally, the application of PHSD to nucleus-nucleus collisions at SPS 
energies has demonstrated a good reproduction of a large set of data [18] improving 
the quality of the description within the Hadron-String Dynamics (HSD) approach 
which lacks explicit partonic degrees of freedom. These general properties of PHSD 
results for idealized systems are well in line with global observations of experiments at 
RHIC energies [6], however, the actual question is about the description of various 
experimental differential observables in particular at top RHIC energies where a 
rather precise experimental control is possible. 

The paper is organized as follows: In Section 2 we briefly review the PHSD approach, 
recall the input parameters and specify the extensions relative to Ref. [18]. Section 
3 is devoted to actual applications for Au + Au collisions at RHIC energies in 
comparison to experimental data from the STAR, PHENIX, BRAHMS and PHOBOS 
collaborations. Section 4 concludes this study with a summary, discussion of open 
problems and an outlook. 



2 The PHSD approach 



The Parton-Hadron-String-Dynamics (PHSD) approach is a microscopic covariant 
transport model that incorporates effective partonic as well as hadronic degrees of 
freedom and involves a dynamical description of the hadronization process from 
partonic to hadronic matter. Whereas the hadronic part is essentially equivalent 
to the conventional HSD approach [32,33] the partonic dynamics is based on the 
Dynamical QuasiParticle Model (DQPM) [15,16,34] which describes QCD properties 
in terms of single-particle Green's functions (in the sense of a two-particle irreducible 
(2PI) approach). In Ref. [18] we have fitted the (essentially three) DQPM parameters 
for the temperature-dependent effective coupling to the lattice QCD results of Ref. 
[35] which lead to a critical temperature T c ps 192 MeV that corresponds to a critical 
energy density of e c ~ 1.25 GeV/fm 3 . These lattice QCD results disagree with the 
1QCD 'data' by the Wuppertal-Budapest group [36] - as pointed out explicitly in the 
summary of Ref. [18] - and the conflict meanwhile has come close to an end [37]. 

In this respect in the present study we refitted the DQPM parameters for the strong 

1 For dissipative hydrodynamics with a finite shear viscosity r\ we refer the reader to Refs. 
[23-31]. 
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coupling in order to reproduce the lattice QCD results from Ref. [38]. The latter 
yield a critical temperature T c ps 160 MeV and a critical energy density e ~ 0.5 
GeV/fm 3 which are significantly lower than those deduced from the results of Ref. 
[35]. Furthermore, the scaled interaction measure (e — 3P)/T 4 (with P denoting the 
pressure) is also significantly lower in the 1QCD calculations from Ref. [38] than in the 
calculations from Cheng et al. [35]. As will be shown below, the readjustment of the 
DQPM parameters to the 1QCD data of Ref. [38] results in a smaller dynamical width 
of the quasiparticles, whereas the dynamical pole masses become slightly larger. For 
the actual procedure of fixing the effective (temperature dependent) strong coupling 
g(T/T c ) we refer the reader to Section 2.1 of Ref. [18]. 

One might worry that the quasiparticle properties - fixed in thermal equilibrium 
- also should be appropriate for out-off equilibrium configurations. This question is 
nontrivial and can only be answered by detailed model investigations e.g. on the basis 
of Kadanoff-Baym equations. We recall that such studies have been summarized in 
Ref. [14] for strongly interacting scalar fields that initially are far off-equilibrium and 
simulate momentum distributions of colliding systems at high relative momentum. 
The results for the effective parameters M and 7, which correspond to the time- 
dependent pole mass and width of the propagator, indicate that the quasiparticle 
properties - except for the very early off-equilibrium configuration - are close to 
the equilibrium mass and width even though the phase-space distribution of the 
particles is far from equilibrium (cf. Figs. 8 to 10 in Ref. [14]). Accordingly, we will 
adopt the equilibrium quasiparticle properties also for phase-space configurations 
out-off equilibrium as appearing in relativistic heavy-ion collisions. The reader has 
to keep in mind that this approximation is far from being arbitrary, however, not 
fully equivalent to the exact solution. 



2.1 Quasiparticle properties and thermodynamics within the DQPM 



The actual gluon mass M g and width 7 9 - employed as input in the further calcula- 
tions - as well as the quark mass M q and width 7^ are depicted in Fig. 1 as a function 
of T/T c . These values for the masses are slightly larger than those presented in Ref. 
[18] and the width 7^ as well as the width 7^ are smaller due to the lower scaled 
interaction measure in Ref. [38] compared to the interaction measure in Ref. [35]. 
This implies that the partons become better 'quasiparticles' since the ratios 7 9 /M g 
and 7 9 /M 9 decrease relative to Ref. [18]. Note that for /j, q — the DQPM gives 

M q = 2 -M g) lq = . (1) 



These variations of the masses with the temperature T that appear drastic in Fig. 1 
become, however, rather smooth if viewed as a function of the scalar parton density 
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Fig. 1. The effective gluon mass M g and witdh j g as function of the scaled temperature 
T/T c (red lines). The blue lines show the corresponding quantities for quarks. 

p s defined (in thermal equilibrium) by 

p,(|) = d J^Jj^ 2 Jp 2 p) M"/t) e( P 2 ) (2) 



/doj f d?p / — 
2^ J (27T? 2 v ^ P<? ^' ^ ~ 

where and rip denote the Bose and Fermi functions, respectively, while p q stands 
for the quark chemical potential. The number of transverse gluonic degrees of freedom 
is d g = 16 while the fermion degrees of freedom amount to d q = dq = 2N c Nf = 18 
in case of three flavors (Nf=3). The function Q(p 2 ) (with p 2 = u 2 — p 2 ) projects on 
time-like four-momenta since only this fraction of the four-momentum distribution 
can be propagated within the light cone. In Eq. (2) the parton spectral functions pj 
are no longer 5— functions in the invariant mass squared but taken as [4] 



P>) -^-((a;-4) 2 + 7l (^ + 4) 2 + 7j) (3) 

separately for quarks and gluons (j = q,q,g). With the convention E 2 (p) = p 2 + 
M 2 — 7J, the parameters M 2 and jj are directly related to the real and imaginary 
parts of the retarded self-energy, e.g. IT,- = M 2 — 2ijjU. The spectral function (3) is 
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Fig. 2. The effective gluon mass M g and width j g as function of the scalar density p s within 
the DQPM. The corresponding mass and width for quarks (for p q = 0) is obtained from 
Eq. (1). Note the logarithmic scale in p s . 

antisymmetric in uj and normalized as 



The dependence of the gluon mass M g and width 7 5 as a function of p s (within the 
DQPM) is displayed in Fig. 2 and demonstrates that the explicit variation with p s is 
rather moderate in view of the logarithmic scale in p s . Note that in transport theory 
the scalar forces on a 'particle' are given by the ratio of the particle mass over its 
energy times the gradient of the scalar mean-field U s (x). The latter gradient is con- 
ventionally written as VU s (x) = dU s /dp s Vp s (x) which demonstrates the separation 
of geometry - expressed by Vp s (x) - from the strength of the force determined by 
dU s /dp s . 

With the quasiparticle properties (or propagators) fixed (cf. Fig. 1 and Eq. (3)) one 
can evaluate the entropy density s(T), the pressure P{T) and energy density e(T) in 
a straight forward manner (cf. Ref. [18]). A direct comparison of the resulting entropy 
density s(T) and energy density e(T) from the DQPM with 1QCD results from Ref. 
[38] is presented in Fig. 3. Both results have been divided by T 3 and T 4 , respectively, 
to demonstrate the scaling with temperature. We briefly note that the agreement 
is sufficiently good. This also holds for the dimensionless 'interaction measure', i.e. 
(e — 3P)/T 4 as demonstrated in Fig. 4. 

We recall that the DQPM uniquely defines a potential energy density 
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V P (T, p q ) = T™(T, p q ) + l£(T, p q ) + T°°(T, p q ) 



(5) 
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Fig. 3. The scaled entropy density s(T)/T 3 (blue line) and scaled energy density e(T)/T 
(red line) from the DQPM in comparison to the 1QCD results from Ref. [38] (full dots and 
triangles) . 
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Fig. 4. The dimensionless 'interaction measure' (e — 3P)/T 4 within the DQPM in compar- 
ison to the 1QCD results from Ref. [38] (full dots). 

where the different contributions correspond to the space-like part of the energy- 
momentum tensor component Tj° of parton j = g,q,q (cf. Section 3 in Ref. [16]). 
As demonstrated in Ref. [16] this quantity is practically independent on the quark 
chemical potential (for moderate p q ) when displayed as a function of the scalar 
density p s instead of T and p q separately. Note that the field quanta involved in (5) 
are virtual and thus correspond to partons exchanged in interaction diagrams. 

A scalar mean- field U s (p s ) for quarks and antiquarks can be defined by the derivative, 

dV p (p s ) 



U s ( Ps ) 



dps 



(6) 
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Fig. 5. The scalar mean field (6) for quarks and antiquarks from the DQPM as a function 
of the scalar parton density p s (2). 

which is evaluated numerically within the DQPM. The actual result for the new 
parameter-set is displayed in Fig. 5 as a function of the parton scalar density p s and 
shows that the scalar mean field is in the order of a few GeV for p s > 10 fm -3 . The 
mean-field (6) is employed in the PHSD transport calculations and determines the 
force on a quasiparticle j, i.e. ~ Mj/EjVU s (x) = Mj/Ej dU s /dp s Vp s (x) where the 
scalar density p s (x) is determined numerically on a space-time grid (see below). 

We point out that in general the quasiparticle masses Mj as well as the widths 7j 
should depend also on the four-momentum q relative to the medium at rest and 
approach the perturbative values at high q 2 . So far, the momentum-dependence of 
the complex self energy cannot reliably be extracted from the 1QCD results in ther- 
modynamic equilibrium which are essentially sensitive to momenta in the order of 
a few times the temperature. This is presently an open issue and will have to be 
re-addressed in future. 



2.2 Hadronization 

The hadronization, i.e. the transition from partonic to hadronic degrees of freedom, 
was presented in Refs. [17,18] and is described below again in more detail for clarity. 
The hadronization is implemented in PHSD by local covariant transition rates as 
introduced first in Ref. [17] e.g. for q+q fusion to a mesonic state m of four-momentum 
p = (lu, p) at space-time point x = (t, x): 




= Tr q Trq 5 4 (p - p g - p q ) 5 4 ( Xq ~^ Xq - x^j 



XW, p q {p q ) Ug p q (Pq) \Vqq\ 2 W m (x q ~ Xq, (j) q ~ Pq)/2) 
xNg(Xq,p q ) Nq(Xq,Pq) 5(flaVOr, Color) . 



(7) 



In Eq. (7) we have introduced the shorthand notation, 



where J2j denotes a summation over discrete quantum numbers (spin, flavor, color); 
Nj(x,p) is the phase-space density of parton j at space-time position x and four- 
momentum p. In Eq. (7) ^(flavor, color) stands symbolically for the conservation of 
flavor quantum numbers as well as color neutrality of the formed hadronic state m 
which can be viewed as a color-dipole or 'pre-hadron'. Furthermore, v q q(p p ) is the 
effective quark- ant iquark interaction from the DQPM (displayed in Fig. 10 of Ref. 
[16]) as a function of the local parton (q + q + g) density p p (or energy density). 
Furthermore, W m (x,p) is the dimensionless phase-space distribution of the formed 
'pre-hadron', i.e. 

W m (&P«) = exp (^) exp (2b\p\ - (M q - M,) 2 /4)) (9) 



with £ = x\ — X2 = x q — Xq and p^ = {p\ — p^jl = (p g — Pg)/2 (which had been 
previously introduced in Eq. (2.14) of Ref. [39]). The width parameter b is fixed by 
(r 2 ) = b = 0.66 fm (in the rest frame) which corresponds to an average rms radius of 
mesons. We note that the expression (9) corresponds to the limit of independent har- 
monic oscillator states and that the final hadron-formation rates are approximately 
independent of the parameter b within reasonable variations. By construction the 
quantity (9) is Lorentz invariant; in the limit of instantaneous 'hadron formation', 
i.e. £° = 0, it provides a Gaussian dropping in the relative distance squared (ri — r2) 2 . 
The four- momentum dependence reads explicitly (except for a factor 1/2) 

{E\ - E 2 f - ( Pl - p 2 ) 2 - (M 1 - M 2 f < (10) 



and leads to a negative argument of the second exponential in (9) favoring the fusion 
of partons with low relative momenta p q — pq = p\ — p 2 - 

Related transition rates (to Eq. (7)) are defined for the fusion of three off-shell quarks 
{qi + Q2 + 93 "H- B) to a color neutral baryonic (B or B) resonances of finite width 
(or strings) fulfilling energy and momentum conservation as well as flavor current 
conservation (cf. Section 2.3 in Ref. [18]). In contrast to the familiar coalescence 
models [40,41] and their recent extensions [42-45] this hadronization scheme solves 
the problem of simultaneously fulfilling all conservation laws and the constraint of 
entropy production. 



9 



2.3 Numerical aspects 



On the hadronic side PHSD includes explicitly the baryon octet and decouplet, the 
0~- and l~-meson nonets as well as selected higher resonances as in HSD [32,33]. 
Hadrons of higher masses (> 1.5 GeV in case of baryons and > 1.3 GeV in case 
of mesons) are treated as 'strings' (color-dipoles) that decay to the known (low- 
mass) hadrons according to the JETSET algorithm [46]. We discard an explicit 
recapitulation of the string formation and decay and refer the reader to the original 
work [46] or Ref. [47]. 

The dynamical evolution of the system is entirely described by the transport dynam- 
ics in PHSD incorporating the off-shell propagation of the partonic quasiparticles 
according to Refs. [11,12,14] as well as the transition to resonant hadronic states (or 
'strings') via Eq. (7). The time integration for the testparticle-equations of motion 
(cf. Refs. [12]) is performed in the same way as in case of hadronic off-shell transport 
where (in view of the presently momentum-independent width 7) the simple relation 
(19) in Ref. [48] is employed. For the collisions of partons two variants are at our 
disposal: i) geometrical collision criteria as employed in standard hadronic trans- 
port, ii) the in-cell method developed in Ref. [49]. The latter can easily be extended 
to describe 2 <H> 3 processes etc. in a covariant way [50]. It is the better choice at 
high particle densities (cf. Ref. [51]) and was actually used in the calculations pre- 
sented below. The hadronization is performed by integrating the rate equations (e.g. 
(7)) in space and time which are discretized on a four- dimensional grid by At and 
AV(t) = Ax(t)Ay(t)Az(t). In beam direction we use an initial grid size Az = l/^ cm 
fm with 7 cm denoting the Lorentz-7 factor in the nucleon-nucleon center-of-mass 
system while in the transverse direction we use Ax = Ay = 1 fm. The grid size 
is increased dynamically during the transport calculation such that all particles are 
included on the actual grid. This practically implies that the grid boundary in beam 
direction approximately moves with the velocity of light. In each time step At and 
cell AV the integrals in (7) and the respective integrals for baryon (antibaryon) for- 
mation are evaluated by a sum over all (time-like) testparticles using (e.g. for the 
quark density) 



3C 



AV J J 2tt q J (2tt) 

AV 



where the sum over J q implies a sum over all testparticles of type q (here quarks) in 
the local volume AV in each parallel run. In Eq. (11) iV denotes the occupation num- 
ber in phase space which in thermal equilibrium is given by Bose- or Fermi-functions, 
respectively. In case of other operators like the scalar density, energy density etc. the 
number 1 in Eq. (11) has to be replaced by ^Jpj/ujj, Uj etc. In order to obtain lower 
numerical fluctuations the integrals are averaged over the parallel runs (typically 50 
at RHIC energies). For each individual testparticle (i.e. x q and p q fixed) the addi- 
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tional integrations in (7) give a probability for a hadronization process to happen; 
the actual event then is selected by Monte Carlo. Energy-momentum conservation 
fixes the four-momentum p of the hadron produced and its space-time position x is 
determined by (7). The final state is either a hadron with flavor content fixed by 
the fusing quarks (and/or antiquarks) or by a string of invariant mass \/s (with the 
same flavor), if y/s is above 1.3 GeV for mesonic or above 1.5 GeV for baryonic quark 
content. 

On the partonic side the following elastic and inelastic interactions are included 
in PHSD qq -B- qq, qq -B- qq, gg B- gg, gg B- g, qq B- g exploiting 'detailed- 
balance' with interaction rates again from the DQPM [16,18]. Numerical tests of the 
parton dynamics with respect to conservation laws, interaction rates in and out-off 
equilibrium in a finite box with periodic boundary conditions have been presented in 
Ref. [52]. For further details we refer the reader to Section 2.2 of [18]. The interactions 
between hadrons are the same as in the HSD transport model. 

2.4 Initial conditions 

The initial conditions for the parton/hadron dynamical system have to be specified 
additionally. In order to describe relativistic heavy-ion reactions we start with two 
nuclei in their 'semi-classical' groundstate, boosted towards each other with a velocity 
f5 (in z-direction), fixed by the bombarding energy. The initial phase-space distribu- 
tions of the projectile and target nuclei are determined in the local Thomas-Fermi 
limit as in the HSD transport approach [53,32,33] or the UrQMD model [54,55]. We 
recall that at relativistic energies the initial interactions of two nucleons are well 
described by the excitation of two color-neutral strings which decay in time to the 
known hadrons (mesons, baryons, antibaryons) [46]. Initial hard processes - i.e. the 
short-range high-momentum transfer reactions that can be well described by pertur- 
bative QCD - are treated in PHSD (as in HSD) via PYTHIA 5.7 [56]. The novel 
element in PHSD (relative to HSD) is the 'string melting concept' as also used in the 
AMPT model [57] in a similar context. However, in PHSD the strings (or possibly 
formed hadrons) are only allowed to 'melt' if the local energy density e(x) (in the 
local rest frame) is above the transition energy density e c . The present DQPM ver- 
sion (fitted to the 1QCD results from Ref. [38]) gives e c ~ 0.5 GeV/fm 3 . The mesonic 
strings then decay to quark-antiquark pairs according to an intrinsic momentum 
distribution, 

F(q) ~ exp(-26 2 q 2 ) , (12) 

in the meson rest-frame (cf. Eq. (7) for the inverse process). The parton final four- 
momenta are selected randomly according to the momentum distribution (12) (with 
b= 0.66 fm), and the parton-energy distribution is fixed by the DQPM at given energy 
density e(p s ) in the local cell with scalar parton density p s . The flavor content of the 
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qq pair is fully determined by the flavor content of the initial string. By construction 
the 'string melting' to massive partons conserves energy and momentum as well as 
the flavor content. In contrast to Ref. [57] the partons are of finite mass - in line with 
their local spectral function - and obtain a random color c = (1,2,3) or (r,b,g) in 
addition. Of course, the color appointment is color neutral, i.e. when selecting a color 
c for the quark randomly the color for the antiquark is fixed by — c. The baryonic 
strings melt analogously into a quark and a diquark while the diquark, furthermore, 
decays to two quarks. 



2. 5 Shadowing 



As well known from deeply-inelastic lepton scattering on nuclei, there is a depletion 
of low-momentum partons in a nucleon embedded in a nucleus compared to the pop- 
ulation in a free nucleon, which leads to a lowering in the production cross section 
of hard probes in proton-nucleus and nucleus-nucleus collisions at high bombarding 
energies compared to the production in the superposition of independent nucleon- 
nucleon collisions. The reasons for depletion, though, are numerous, and models of 
shadowing vary accordingly. There is, therefore, a considerable (about a factor of 3) 
uncertainty in the amount of shadowing predicted at RHIC and especially at LHC 
[58-65] . In the analysis of the d + Au data at y / J/v r /v = 200 GeV, in which the max- 
imum estimate for the effect of the shadowing was made [59,65,66], the shadowing 
lead to a ~ 10% reduction while an anti-shadowing closer to target rapidities was 
observed. More recent estimates for the shadowing at RHIC energies have been pre- 
sented in Refs. [63,64] and imply a lower amount of shadowing for hard probes. In 
the present PHSD calculations we employ the results of Ref. [58] which lead to a 
suppression of charmonia at forward rapidities in d+Au reactions due to shadowing 
by less than 5%. 

The actual implementation of shadowing in PHSD is done in a practical way by 
parametrizing the suppression (or enhancement) factors from Ref. [58] for heavy 
nuclei at RHIC energies as a function of the Bjorken variable x ~ 2p /y/s and the 
mass number A, i.e. 



R A (x) = 1 + 1.19 ln 1/6 [x 3 - 1.5(x + x L )x 2 + 3x x L x] 
1.08(^/3 _ 



(13) 



ln(A + 1) 



exp(-a; /x ), 



with the parameters x = 0.1, x L = 0.7 and a a = 0.1 (A 1 ^ — 1). Since Eq. (13) 
determines only the average nuclear effect one has to specify the impact parameter 
dependence of the shadowing function Ra(x) on the actual position of the colliding 
nucleons. As in Ref. [58] we assume that the shadowing parameter a a is proportional 
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to the longitudinal thickness of the nucleus at impact parameter r and adopt 



where r is the transverse distance of the interacting nucleon from its nucleus center 
and Ra is the radius of its nucleus. In this manner an approximate implementa- 
tion of initial state shadowing is achieved which is sufficient for the purposes of the 
present investigations. Since the shadowing effects are only on the level of a few per- 
cent for the observables addressed here we discard a further discussion and shift its 
representation to a forthcoming study of hard probes. 

Let's summarize the modifications and extensions of the PHSD approach relative to 
the version presented before: except for shadowing and the new partonic equation- 
of-state from Ref. [38], we adopt here essentially the same approach as described in 
detail in Ref. [18], which we use now to investigate heavy- ion reactions at the top 
RHIC energy while in Ref. [18] the SPS energy regime was studied. 



3 Application to Au + Au collisions at y/s = 200 GeV 

In this Section we employ the PHSD approach - described briefly in Section 2 and 
in more detail in Ref. [18] - to nucleus-nucleus collisions at ultra-relativistic energies, 
i.e. in particular at the top RHIC energy. Note that at RHIC or more specifically 
LHC energies other initial conditions (e.g. a color-glass condensate [67]) might be 
necessary. In the present work we discard such alternative initial conditions and 
explore to what extent the present initial conditions are compatible with differential 
measurements by the various collaborations at RHIC. 

3. 1 Parton dynamics at RHIC energies 

We start with a consideration of energy partitions in order to map out the fraction of 
partonic energy in time for relativistic nucleus-nucleus collisions. In Ref. [18] we have 
found that even in central collisions of Pb +Pb at 158 A GeV only a limited fraction 
of degrees-of-freedom can be attributed to a partonic phase due to a significant 
hadronic corona [68,69] both in coordinate space as well as for large rapidities. In 
order to address the experimental observations at RHIC we will focus here on the 
total energy at midrapidity, i.e. for \y\ < 1. We note in passing that the total energy 
- integrated over all rapidities - is conserved throughout the reaction within better 
than 1% (cf. also Ref. [18]). 

In Fig. 6 we show the energy balance for a central (impact parameter 6=1 fm) reaction 
of Au+Au at y/s = 200 GeV, i.e. at the top RHIC energy including partonic, mesonic 
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Fig. 6. The total energy Eta (upper solid black line) for central (6=1 fm) collisions of 
Au+Au at y/s = 200 GeV within the rapidity window \y\ < 1. The dashed (blue) line 
shows the energy contribution from partons while the dot-dashed (red) line displays the 
energy contribution from mesons (including 'unformed mesons' in strings). The long-dashed 
(green) line is the contribution of baryons (and antibaryons). 

and baryonic degrees of freedom in the rapidity window \y\ < 1. The total energy 
E tot within this rapidity interval (upper line) - which at t = is zero in the cms of the 
colliding nuclei - shows a rapid increase in time for £ ~ 1.6 fm/c which corresponds to 
the contact time of the colliding heavy ions. At t ~ 2.5 fm/c about 85% of the energy 
(at midrapidity) is carried by the partonic degrees-of-freedom which are converted 
with increasing time to mesons and baryons (or antibaryons) essentially within 6- 
8 fm/c. The total energy within this rapidity window is not conserved since by 
elastic and inelastic reactions the reaction products may leave (or enter) the rapidity 
window. Note that even in central collisions (at midrapidity) not all the energy is 
converted to a partonic phase and a hadronic (or rather string-like) corona [68,69] 
survives in the surface area of the collision zone. 



We note in passing that a qualitatively similar picture is obtained when plotting the 
parton, meson and baryon (+ antibaryon) numbers for \y\ < 1. An essential point 
here is that the number of final hadronic states is larger than the maximal number of 
partons, i.e. there is a production of entropy in the hadronization process as pointed 
out before in Refs. [17,18] and thus no violation of the second law of thermodynamics 
in PHSD! 

In order to shed more light on the hadronization process in PHSD we display in Fig. 7 
the invariant mass distribution of qq pairs (solid line) as well as qqq (and qqq) triples 
(dashed line) that lead to the formation of final hadronic states. The reaction is again 
Au + Au at y/s = 200 GeV at impact parameter b— 1 fm. In fact, the distribution 
for the formation of baryon (antibaryon) states starts above the nucleon mass and 
extends to high invariant mass covering the nucleon resonance mass region as well 
as the high-mass continuum (which is treated by the decay of strings within the 
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Fig. 7. The invariant mass distribution for fusing qq pairs (solid blue line) as well as qqq 
(and qqq) triples (dashed red line) that lead to the formation of final hadronic states for 
Au+Au at \/s = 200 GeV (b=l fm). The vertical arrows denote the thresholds for mesonic 
(1.3 GeV) and baryonic strings (1.5 GeV). 

JETSET model [46,47] for M > 1.5 GeV). On the 'pre-mesonic' side the invariant- 
mass distribution starts above the pion mass and extends up to continuum states 
of high invariant mass (described again in terms of string excitations for M > 1.3 
GeV). The low- mass sector is dominated by a, p, a 1; oj or K*, K* transitions etc. As 
mentioned before the excited 'pre-hadronic' states decay to two or more 'pseudoscalar 
octet' mesons such that the number of final hadrons is larger than the initial number 
of fusing partons. 

The individual reaction rates in the partonic phase are of further interest in central 
Au+Au collisions at the top RHIC energy. To this aim we display in Fig. 8 the 
interaction rates for the channels qq — >■ g (black dotted line), g — > qq (dashed red 
line) and elastic parton scattering (dot-dashed green line). Except for the very early 
phase, where the elastic scattering channels dominate, all interaction rates are of 
comparable size and decrease rapidly in time. For comparison we also show the 
hadronization rate by the solid blue line which has a maximum at about 6 fm/c 
and is sizeably larger than the other interaction rates for t > 7 fm/c. Note that 
in this representation we have considered all rapidities; accordingly interaction and 
hadronization processes keep on going at forward and backward rapidities also for 
times larger than 50 fm/c. 

Some information on the time evolution of the quark and gluon mass functions is dis- 
played in Fig. 9 which shows the number of 'particles' as a function of invariant mass 
M and time t at midrapidity < 1). Note that by integration over M one obtains 
the number of quarks (+ antiquarks) N q (t) and gluons N g (t) in the rapidity interval 
\y\ < 1 while dividing by N q (t) and N g (t), respectively, a rough estimate for the par- 
ticle spectral functions is obtained. Note that the mass distributions displayed here 
are the product of the spectral functions and the occupation numbers in a restricted 
phase space. Due to a moderate variation of the partons pole mass and width with 
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Fig. 8. The interaction rates (integrated over rapidity) for the channels qq — >■ g (black 
dot-dashed line), g — )• qq (dashed red line) and elastic parton scattering (dashed green 
line). The solid blue line displays the hadronization rate. The reaction is Au+Au at yfs = 
200 GeV (6=1 fm). The rateinclude all rapidities. 



the scalar density p s (cf. Fig. 2) the shapes of the partonic mass distributions do 
not change very much in time. The average quark mass is about 0.5 GeV while the 
average gluon mass is only slightly less than 1 GeV. Note, however, that the width 
of the mass function - which reflects the actual interaction rate per parton - remains 
significant for all times up to hadronization. Furthermore, the average parton width 
as a function of time cannot directly be related to Figs. 1 or 2 since at a given time 
t the partonic mass distribution relates to different scalar densities in the course of 
the partonic evolution. 




Fig. 9. The time-dependent mass distributions for quarks (+ antiquarks) (l.h.s.) and gluons 
(r.h.s.) for a central Au+Au collision at y/s = 200 GeV and 6=1 fm at midrapidity (|y| < 1). 
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3.2 Particle spectra in comparison to experiment 



Apart from the more general considerations in the previous Subsection, it is of inter- 
est how the PHSD approach compares to the HSD model (without explicit interacting 
partonic degrees of freedom) as well as to experimental data from the RHIC collabo- 
rations. We start with rapidity spectra from PHSD (solid red lines) for charged pions 
and kaons in 5% central Au+Au collisions at = 200 GeV which are compared in 
Fig. 10 to the data from the RHIC Collaborations [70-72] as well as to results from 
HSD (dashed blue lines). We find the rapidity distributions of the charged mesons 
to be slightly narrower than those from HSD and actually closer to the experimental 
data. Also note that there is slightly more production of K ± mesons in PHSD than 
in HSD while the number of charged pions is slightly lower. The actual deviations 
between the PHSD and HSD spectra are not dramatic but more clearly visible than 
at SPS energies (cf. Ref. [18]). Nevertheless, it becomes clear from Fig. 10 that the 
energy transfer in the nucleus-nucleus collision from initial nucleons to produced 
hadrons - reflected dominantly in the light meson spectra - is rather accurately de- 
scribed by PHSD. Fig. 10 also demonstrates that the longitudinal motion is well 
understood within the PHSD approach. 

Independent information on the active degrees of freedom is provided by transverse 
mass spectra of the hadrons especially in central collisions. The actual results for 
RHIC energies are displayed in Fig. 11 where we show the transverse mass spectra 
of 7t _ , K + and K~ mesons for 5% central Au+Au collisions at y/s = 200 GeV in 
comparison to the data of the RHIC Collaborations [70-72]. Here the slope of the 
7r~ spectra is slightly enhanced in PHSD (solid red lines) relative to HSD (dashed 
blue lines) which demonstrates that the pion transverse mass spectra also show 
some sensitivity to the partonic phase (contrary to the SPS energy regime). The 
K ± transverse mass spectra are substantially hardened with respect to the HSD 
calculations - i.e. PHSD is more in line with the data - and thus suggest that partonic 
effects are better visible in the strangeness degrees-of-freedom. The hardening of 
the kaon spectra can be traced back to parton-parton scattering as well as a larger 
collective acceleration of the partons in the transverse direction due to the presence of 
the repulsive mean-field for the partons (cf. Fig. 5). The enhancement of the spectral 
slopes for kaons and antikaons in PHSD due to collective partonic flow shows up 
much clearer for the kaons due to their significantly larger mass (relative to pions). 

The latter considerations also become transparent when comparing the transverse 
mass spectra for protons at midrapidity from HSD and PHSD to the data from 
the PHENIX Collaboration [70] in Fig. 12. Here the HSD spectra (dashed blue line) 
severely underestimate the slope of the data from Ref. [70] whereas the PHSD spectra 
(solid red line) are fairly in line. These differences are so dramatic because in HSD 
the protons at midrapidity dominantly stem from initial string decays and are not 
allowed to rescatter during their formation time of 7lTo where 7^ denotes the Lorentz 
factor and r = 0.8 fm/c is the default formation time for hadrons in HSD [32,33]. On 
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Fig. 10. The rapidity distribution of tt + (upper part, l.h.s.), K + (lower part, l.h.s.), 7r~ 
(upper part, r.h.s.) and K~ (lower part, r.h.s.) for 5% central Au+Au collisions at yfs = 
200 GeV from PHSD (solid red lines) in comparison to the distribution from HSD (dashed 
blue lines) and the experimental data from the RHIC Collaborations [70-72]. 



the other hand, in PHSD the dominant source of protons (and also the other baryons 
and antibaryons) at midrapidity is the fusion of three quarks from the partonic phase. 
Since the partonic degrees of freedom interact strongly and are accelerated in the 
expansion phase due to the repulsive mean field the protons pick up the momenta 
from the fusing partons and thus show a sizeable harder slope in the transverse mass 
spectrum. 



In summarizing this Subsection we point out that the partonic phase in PHSD at the 
top RHIC energy leads to a narrowing of the longitudinal momentum distribution, 
a reduction of pion production and slight enhancement of kaon production and to 
a hardening of their transverse mass spectra relative to HSD (closer to the data). 
These effects are clearly visible especially in the transverse degrees-of-freedom and 
are identical to those of Ref. [18] at SPS energies, however, more pronounced due to 
the larger space-time region of the partonic phase. 
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Fig. 11. The 7r~, and K~ transverse mass spectra for 5% central Au+Au collisions at 
yfs = 200 GeV from PHSD (solid red lines) in comparison to the distributions from HSD 
(dashed blue lines) and the experimental data from the BRAHMS, PHENIX and STAR 
Collaborations [70-72] at midrapidity. 
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Fig. 12. The proton transverse mass spectra for 5% central Au+Au collisions at yfs = 200 
GeV from PHSD (solid red line) in comparison to the distributions from HSD (dashed blue 
line) and the experimental data from the PHENIX Collaboration [70] at midrapidity. 

3.3 Elliptic flow 

Of additional interest are the collective properties of the strongly interacting system 
which are explored experimentally via the elliptic flow 

V 2 (PT, V) = ({pi - Pi)/ (Pi + Vl)) \pt, V (15) 
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of hadrons as a function of centrality, rapidity y, transverse momentum px or trans- 
verse kinetic energy per participating quarks and antiquarks. We note that the re- 
action plane in PHSD is given by the x — z plane with the z-axis in beam direction. 

We start in Fig. 13 with the elliptic flow v 2 (for Au+Au collisions at the top RHIC 
energy) as a function of the centrality of the reaction measured by the number of 
participating nucleons N part . The solid (red) line stands for the results from PHSD 
which is compared to the data for charged particles from the PHOBOS Collaboration 
[73] . The dashed blue line refers to the corresponding results for v 2 from HSD (taken 
from Ref. [75]). The momentum integrated results in the pseudo-rapidity window 
\r)\ < 1 from PHSD compare well to the data from Ref. [73] whereas the HSD results 
clearly underestimate the elliptic flow as pointed out before in Ref. [75]. The relative 
enhancement of v 2 in PHSD with respect to HSD can be traced back to the high 
interaction rate in the partonic phase and to the repulsive mean field for partons (cf. 
Fig. 5). We note in passing that PHSD calculations without mean fields only give a 
tiny enhancement for the elliptic flow relative to HSD. 
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Fig. 13. The elliptic flow v 2 for Au+Au collisions at the top RHIC energy yfs = 200 GeV as 
a function of the centrality measured by the number of participating nucleons N par t- The 
solid (red) line stands for the results from PHSD whereas the dashed (blue) line represents 
the results from HSD (from Ref. [75]). The data are taken from the PHOBOS Collaboration 
[73] and correspond to momentum integrated events in the pseudo-rapidity window |r/| < 1 
for charged particles. The shaded band signals the statistical uncertainties of the PHSD 
calculations. 

Fig. 14 shows the final hadron v 2 versus the transverse momentum pt for different 
particle species in comparison to the data from the STAR [74] and PHENIX Collab- 
orations [8]. We observe a mass separation in px as well as a separation in mesons 
and baryons for p? > 2 GeV roughly in line with data. The elliptic flow of mesons is 
slightly underestimated for p? > 2 GeV in PHSD which is opposite to ideal hydro- 
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Fig. 14. The hadron elliptic flow V2 for inclusive Au+Au collisions as a function of the 
transverse momentum px (in GeV) for different hadrons in comparison to the data from 
the STAR [74] and PHENIX Collaborations [8] within the same rapidity cuts. 

dynamics which overestimates v<i at high transverse momenta. On the other hand, 
the proton (and antiproton) elliptic flow is slightly overestimated at low transverse 
momenta pr < 1-5 GeV. 

A further test of the PHSD hadronization approach is provided by the 'constituent 
quark number scaling' of the elliptic flow V2 which has been observed experimentally 
in central Au+Au collisions at RHIC [6,8,9]. In this respect we plot V2/n q versus the 
transverse kinetic energy per constituent parton, 



KE 7 



ttit — rn 
n„ 



(16) 



with rriT and m denoting the transverse mass and actual hadron mass, respectively. 
For mesons we have n q = 2 and for baryons/antibaryons n q = 3. The results for the 
scaled elliptic flow are shown in Fig. 15 in comparison to the data from the STAR [74] 
and PHENIX Collaborations [8] for different hadrons and suggest an approximate 
scaling. For KEt > 0.5 GeV there is a tendency to underestimate the experimental 
measurements for A, S,A, S baryons which we attribute to an underestimation of 
interaction terms in PHSD for high momentum hadrons. In this respect we recall 
that the momentum independence of the quasiparticle width 7 and mass M (cf. 
Subsection 2.1) is presently a rough approximation and has to be refined. Due to the 
limited statistics especially in the baryonic sector with increasing pt this issue will 
have to be re-addressed with high statistics in future. 
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Fig. 15. The elliptic flow V2 - scaled by the number of constituent quarks n q - versus the 
transverse kinetic energy (16) divided by n q for different hadron species in comparison to 
the data from the STAR [74] and PHENIX Collaborations [8]. 

4 Summary and outlook 

In summary, relativistic collisions of Au+Au at top RHIC energies have been studied 
within the PHSD approach which includes explicit partonic degrees of freedom as 
well as dynamical local transition rates from partons to hadrons (cf. Eq. (7)). The 
partonic equation-of-state employed has been adopted from the 1QCD calculations of 
the Wuppertal/Budapest group in thermodynamic equilibrium [38] and incorporated 
in the PHSD approach by means of the Dynamical QuasiParticle Model (DQPM) 
[16]. The transition from partonic to hadronic degrees of freedom is described by 
covariant transition rates for the fusion of quark- ant iquark pairs or three quarks 
(antiquarks), respectively, obeying flavor current-conservation, color neutrality as 
well as energy-momentum conservation. 

Our dynamical studies for heavy-ion collisions at relativistic collider energies have 
been compared to earlier results from the Hadron-String Dynamics (HSD) approach - 
incorporating no explicit interacting partonic phase - as well as to experimental data 
from the STAR, PHENIX, BRAHMS and PHOBOS collaborations for Au+Au col- 
lisions at the top RHIC energy of y/s = 200 GeV. We find a reasonable reproduction 
of hadron rapidity distributions and transverse mass spectra and also an acceptable 
description of the elliptic flow of charged hadrons as a function of the centrality of 
the reaction and the transverse momentum p?- This result is quite remarkable since 
the additional interactions of partonic nature (relative to HSD) are essentially deter- 
mined by the DQPM which itself has been fixed by 1QCD 'data'. Furthermore, an 
approximate 'quark-number scaling' of the elliptic flow v 2 is observed in the PHSD re- 
sults, too, while the HSD calculations underestimate the elliptic flow observables. We 
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mention that no fitting has been addressed in our study and that the discrepancies 
with respect to the data of the RHIC collaborations might be attributed to the crude 
approximations with respect to the quasiparticle mass and width, i.e. in particular 
by adopting momentum-independent quantities. This definitely needs improvement 
and more detailed investigations in future. 

Since the bulk dynamics of nucleus-nucleus reactions at RHIC energies appear to 
be reasonably described by the PHSD approach future studies will concentrate on 
leptonic and photonic probes as well as charm and high pt degrees of freedom where 
previous studies within HSD showed sizeable discrepancies with respect to the data 
taken at the SPS or particularly at RHIC energies [76-78] and the necessity for 
partonic degrees of freedom had been pointed out. In fact, the various leading chan- 
nels in the partonic phase for dilepton production have already been calculated on 
the basis of the DQPM propagators in Ref. [79] and been implemented in PHSD. 
Preliminary results have been presented in Ref. [80] and appear encouraging. 

Furthermore, an application of PHSD at LHC energies will be mandatory in order 
to map out the parton dynamics at even higher energy densities. Since the DQPM 
predicts rather moderate changes of the quasiparticle properties when increasing the 
temperature from ~ 2T C to about 4T C only minor changes in the collective properties 
are expected within PHSD when comparing nucleus-nucleus collisions at LHC rela- 
tive to RHIC energies. However, this expectation will have to be studied in detail in 
comparison to the data that have become available recently [81]. 
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